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DESCRIPTION  AND  PURPOSES 


The  level  probabilities  of  order  restricted  inference  are  fundamental 
to  that  theory;  their  values  are  the  probabilities  that  the  order  restricted 
maximum  likelihood  estimates  of  normal  means  assume  specified  numbers  of 
distinct  values,  called  levels.  Those  probabilities  are  computed  under 
the  assumptions  that  the  population  means  are  equal  and  that  the  sampling 
from  the  various  populations  is  done  independently.  The  level  probabili¬ 
ties  depend  upon  a  vector  of  weights  (usually  the  various  sample  sizes) 
and  the  use  of  much  of  the  theory  of  order  restricted  inference  has  been 

limited  by  the  fact  that  those  level  probabilities  can  be  virtually  impos- 

*  \ 

sible  to  compute  if  the  weights  are  not  all  equal. 


Bohrer  and  Chow  (1978;  Algorithm  AS  122)  give  a  program  for  computing 
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these  level  probabilities  when  the  number,  K,  of  populations  is  no  more 
than  10.  Their  program  uses  an  algorithm  for  computing  orthant  probabili¬ 
ties  which  is  due  to  Milton  (1972).  The  time  needed  to  use  Milton's 
algorithm  increases  exponentially  in  K  and  can  require  several  minutes 
or  more  of  computation  when  K  ^  6  (cf.  Bohrer  and  Chow). 

Cran  (1981;  Algorithm  AS  158 )  gives  a  program  for  computing  these  level 
probabilities  when  the  number  K  does  not  exceed  6.  For  K  =  5  it  uses 
an  approximation  due  to  Plackett  (195*0  and  for  K  =  6  an  approximation 
due  to  Childs  (1967)  is  used. 

Robertson  and  Wright  (1982)  develop  an  approximation  which  is  based 
upon  an  idea  of  Chase  ( 197*0  and  uses  the  pattern  of  large  and  small  weights 
We  refer  the  interested  reader  to  Robertson  and  Wright  (1982)  for  an  evalu¬ 
ation  of  the  quality  of  this  approximation.  The  Fortran  program  given 
below  uses  this  approximation  for  the  values  of  the  level  probabilities 
for  K  such  that  Ki  6.  For  K  £  20  and  equal  weights  or  for  general 
weights  and  K  s  5  the  program  is  identical  to  Cran's  (1981;  Algorithm 
AS  158). 
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STRUCTURE 

SUBROUTINE  PROBS  (K,  W,  P,  IFAULT) 


K 

Integer 

input : 

W 

Real  Array  (K) 

input : 

P 

Real  Array  (K) 

output : 

IFAULT 

Integer 

output : 

Auxiliary  Algorithms 
FUNCTION  PRl(l,J,W) 

(Algorithm  AS  158. l) 

FUNCTION  F1(V1,V2,V3) 

(Algorithm  AS  158.2) 

FUNCTION  FACT (M, IFAULT) 

SUBROUTINE  CHASE (K,CH, Pi) 

Chase's  approximations  for  a  given 
SUBROUTINE  PAPRX(K,W,PA,CH,Pl) 

K 
W 
PA 
CH 
PI 


the  number  of  weights 
the  original  weights 
the  computed  probabilities 
a  fault  indicator,  equal  to 

1  if  at  least  one  weight  is  not 
positive 

2  if  K  <  2  or  K  >  20 

3  if  an  error  occurred  in  function 


FACT 

0  otherwise 

computes  explicitly  the  probabilities  for  K  s  5. 

computes  the  correlation  p  = 

computes  n  factorial.  (Algorithm  AS  158.5) 
computes  the  equal  weight  probabilities  and 
K. 

computes  the  approximate  probabilities: 
input:  the  number  of  weights 
input:  the  original  weights 
output:  the  approximate  probabilities 
input:  Chase's  approximate  probabilities 
i..put :  equal  weight  probabilities 


Integer 
Real  Array (K) 
Read  Array (K) 
Real  Array(K,K) 
Real  Array(K,K) 


f  VI  »V3  V 
\  ( V1+V2 ) ( V2+V3 ) / 


Restrictions 


The  weight  array  can  have  no  more  than  20  elements,  so  K  £  20.  In  addi¬ 
tion,  all  weights  must  be  positive. 
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DIMENSION  P  *2  0  >  *  W  (  2  0  ) 


DIMf  NSICN  P  *2  0)  ,  W*20» 

DRIVER  PROGRAM  FOR  SU3R0UTINF  PRCBS. 

READ  <  c. «  10  )  * 

FORMAT (  T  2  ) 

RFAu*  20  )  t  W(l)»  I  =  1*  <  ) 

FORMAT  *  10F».4  ) 

FCHC  CHCFK  THE  INPUT  VALUES. 

UR  I T  r  (  f>,  3  0  > 

FORMAT (  1 H 1  ) 

UR  ITT*  t»  *♦  0  >  K 

FORMAT  (  4  X  *  1  OHTHfP  F  ARE.  «  12*  SHUflGHTS.  ) 

WRITE  (  <jC  > 

FORMAT*  4  X  *  15HTHE  WEIGHTS  ARF  ) 

WRITE  f  t*  tt  >  (  W*D*  I  =  1*  K  ) 

FORMAT*  4  X  *  1  0  F  1  ?  *  '7  ) 

CALL  P*\OHS<  K  *  «*  P*  IF  AULT  ) 

IF  <  I E  AULT  .EG.  1  >  GOTO  ICO 

IF  <  IK  AULT  .EQ.  2  )  GOTO  120 

IF  <  I E  A'JLT  .Eli.  3  >  GOTO  140 

OUTMJT  TKE  COMPUTED  °ROB A3I L IT IES. 

WRITE*  (.«  73  > 

FORMAT*  //4X«  *3HTHF  COMPUTED  PROFAE'IL  ITI  ES  ARE*  /  ) 

DO  oo  L  =  1  *  K 

WPITF*  b *  H  C  )  L*  K  P*L> 

FORMAT*  4  X  »  2HP ( «  12.  1H,  «  12*  4H)  =  «  F10.7  ) 

CONTINUE 

STOP 

WRITE*  F.  110  > 

FORMAT*  4  X  *  35MAT  LEAST  ONE  WEIGHT  IS  NOT  POSITIVE  ) 
STOP 

WRITE*  4*  13C  > 

FORMAT*  4  X  *  17HK  IS  OUT  OF  RANGE  ) 

STOP 

WRITE*  lbO  ) 

FORMAT*  4  X »  34H  AN  ERROR  OCCURRED  IN  FUNCTION  FACT  ) 

STOP 

END 


„  i  to  JiiiU  does  not 

narmit  hilly  legible  leptoduction 


SUBROUTINE  FROOSC  <»  Wt  Pf  IFA'JLT  > 


SUFFCUTTNE  CROP  SC  K  y  W,  P»  TF  AULT  ) 

CALCULATION,  OF  THE  PROP  AGIL  TT  IF  S  PCLtK)  FOR 
THL  C  A  f  <"  Oc  SIMPLE.  ORDER. 

DIMENSION  i»  <  20  >t  P(  20  >t  GC  20 1  20.).  OH  (  2 
DIMENSION  ’■'1C  ?•}.*  20  >♦  PA<  20  .  20  > 

DATA  Cl/  1. OF -6  / 

CHECK  T  H  A  T  "TIGHTS  ARE  POSITJV17. 

IFAl'LT  -  G 
DO  1  I  r  1,  K 

IF  (  U(  I  )  .LL.  C.n  )  GO  TO  Ifil 
CONTINUE 


CHECK  THAT  K  .OC .  *  AND  .LE.  20 


IF 

<  K  .LT. 

5  .  OP. 

K  . 

GT.  20 

) 

10 

T  ^  1? 

2 

UU 

=  U  <  l 

) 

DO 

2  I  = 

2 1 

K 

IF 

(  A  ‘  ~C 

* 

L  -  W  ( 

T 

) 

)  .GT. 

Cl 

) 

GC  TC 

7 

CONT IN.T 

EGO  TL 

‘ft 

eights 

QC 

It  '  ) 

- 

l . :  / 

•? 

.0 

QC 

2 1  3  > 

r 

O.G 

OC 

?t  *■  > 

- 

1.0  / 

f 

.c 

IF 

(  K  .E 

r> 

■  t  • 

3  )  GO 

TO 

H 

DO 

4  J  - 

4* 

K 

AJ 

=  J 

A  1 

=■  1.0 

/ 

AJ 

A? 

=  (  AJ 

- 

1.0  ) 

* 

A  1 

QC 

It  J  ) 

= 

A  1 

Jl 

-  J  - 

1 

no 

3  L  r 

2, 

Jl 

Ll 

=  L  - 

1 

GC 

Lt  J  ) 

- 

A 1  *  G 

( 

Ll 

t  Jl  ) 

♦ 

A  2 

*  QC 

Lt 

CONT  IN.j' 

0(  J.  V  >  -  1.0  /  c  ACT (  J,  I F  AULT  > 

IF  I  TFAULT  .NE.  C  )  RCToRN 
CONT  IN’.''. 

CONT  INI'" 

DO  i>  J  -  1  «  K 
PC  J  )  “  GC  Jt  K  > 

CONT  TNT 
RETURN 

IJNEU'AL  .EIGHTS  -  CHECK  THA  >  , LE.  (■ 


IF  <  K  .GT.  b  )  CO  TO  II 
K2  =  K  -  2 

GO  TC  <  8 «  rt  1C  ).  K? 

PC  i  )  =  PK1C  Ir  3*  W  ) 
PC  2  >  "  0. ‘ 

PC  3  )  -  0.c  -  (’C  1  ) 
RETURN 

PC  1  )  -  p»  1C  It  A.  J  ) 
PC  A  )  *  PR.  1C  'It  4.  U  ) 


SUFfcfiUriNr  r  fr  o  ' c  <  K,  W.  F*  IF  AULT  ) 


pc  2  >  -  o . u  -  p(  4  ) 

p  <  ’  )  -  0 .  c  -  •  (  l  ) 

RETURN 

PC  r>  )  =  PR  1  (  **  5*  W  ) 

P(  4  )  =  PU<  4*  b*  4  > 

P<2)=0.C-PC4> 

P  <  1  >  =  PR  K  1  .  5.  V  ) 

p  <  3  )  =  o .  r.  -  p(i)-pcr) 

RETURN 

CALL  CMSFC  K  t  nHi  pl  ) 

CALL  PA»RX<  K*  M  FA*  CS  PI  > 

no  1  ?  .1  =  1  *  K 

PC  J  )  -  P A (  J ,  K  ) 

CON  T  IM  r 

RCTURN 

IFAULT  -  1 

RETURN 

IFAULT  r  2 

RETURN 

END 

FUNr  TIC.  FF1C  T  «  J,  J  ) 

ALGORITHM  A S  15''.  1  ♦  APPL  STAT  (  1*Tpj  ),  VOL  30,  N 

EXPLICIT  CALCULATION  OF  PRO1  APIL  TTI 1  ?  FOR  K  .LE. 

A  LS  C  CALLED  LiY  FUNCTION  F  ? 


DIMENSION  VC  IT  ) 

DATA  oil/  C  .31 £3  ODES*  / 
ir  (  J  .N't.  3  J  GO  TO  4(1 

C  =  0.!  *  ->11  *  F 1  (  VC  1  ),  '*'(  2  >♦  VC  '  >  ) 

IF  <  I  .EG.  3  )  t»0  TO  30 

PR  J  =  ii  .‘5  b  -  C 

RETURN 

PR  1  -  :.?5  ♦  c 

RETURN 

W  1  =  W  <  1  > 

U2  =  V<  2  ) 

W3  =  3  ) 

W  4  =  *'  (  4  ) 

U1C  =  >. 1  ♦  U? 

U23  =  .0  ♦  l>3 

W34  =  \i  *  \  4 

Slit  -  fl(  U 1 ,  I*  0  *  W 3  ) 

S23  =  c:<  i»r,  v»,  w4  > 

IF  (  J  .E  Q .  b  )  GO  TO  bG 

IF  (  I  ,FS.  4  )  GO  TO  41 

Cl  -  0..’b  *  PIT  * 

CMC.1,W?,U34>  ♦  F1CJI,W?3,V4  )  ♦  F 1 ( .1 ? , » 3 , W 4  )  > 

pr  l  =  ;:  .125  -  r  i 
RETURN 

C2  =  O.^  *  Pll  •  <  SI*  ♦  ::C3  ) 

PR  1  =  0 .125  ♦  C? 

RETURN  c 

V 5  -  VC  5  ) 

W4b  =  Vi  4  ♦  !»b 
Vl?3  =  VI 2  ■*  W3 

W234  =  *23  ♦  W  4 

W34b  =  ^34  •*  V5 


SUC  R  P'JT  IMP  ^0"S(  J,  p,  IFAt.'LT  > 
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S  3R  =  r  i  (  .  y  h  ,  W‘D  ) 

IF  (  1  .EC.  4  )  GC  Tf)  r>? 

C5  =  0 .  05  2 <'  ♦  0.1?r  *  PTI  *  <  M2  ♦  S2  3  ♦  S3H  >  ♦ 

*  0.25  *  pi  i  .  n?  *  s’ a 
IF  (  T  .EO.  1  >  GO  TO  bl 
PP1  -  Cb 

IF  <  P-1  .LT.  ‘>.0  )  PR1  =  0.0 
RETURN 

S  1  1  .*•  p  M<  «I  •  ,!?♦  W3*5  ) 

S151  =  rl(  vi,  ,234.  W‘.  > 

S3 11  =  1(  Wl23.  UR  •  W1'  > 

C3  =  0.M5  ♦  0.125  *  CII  *  (Sll?  ♦  S13J  ♦  SMI  +  Fl<Wl,  W2Z*  »’A5) 

*  ♦  F  i « y  1  i;  ,  y.M  W'tb)  *  Ficwl2*  W3H»  wr>  -  M2  -  s 2 3  -  S3**  ) 

*  -  0.25  *  PII  *  <MI  *  <  SI  2  *  S  3 1  1  ♦  >?3  «  SIM  ♦  S34  *  S113  ) 

PR1  =  0.5  -  03  -  CP 

RETURN 

C?  =  0.125  *  PII  *  ( SI  2  ♦  S3  A  ■*•  Fl(';i,  W2  ♦  W34)  ♦  ri(WI2,  V3t  *4) 

*  ♦  F 1  ( y  1 ,  w?3%  v*A)  ♦  F1<W?*  W?  *  Who)  ♦  F 1  (  W2 i*  WH,  Wt.)  ♦ 

*  F 1  ( W  2  «  bM(  W5  )  ) 

PR  1  =  (!  .2 5  ♦  C2 
RETURN 

FND 

FUNC  ri C\  F 1  (  VI  ,  V?,  V3  ) 

A  L  f>  0  p  I  T  H  AS  lbfi.2  5PPL  ST  AT  <1‘?«1>  VOL  30  t  ?!0  l 

RH2  =  -SfJRTC  VI  *  V3  /  <  (  VI  ♦  V2  )  *  (  V?  ♦  VT  >  )  ) 

FI  =  A5nt  PMO  ) 

RETURN 

END 

FUNCTION  FAf.TC  IFAULT  ) 

ALGCMTHP  AS  lr.«.5  4PPL  ST  AT  C  l^fll  )  VOL  30*  NO  1 

CALCULATION  OF  V  FACTORIAL 

DATA  M A  v M /  bO  / 

IFAULT  =  * 

I F  (  M  .LT.  C  .hr.  y  .  G  T.  'MX*  )  R  r  TUR  f 
IFAIJLT  -  0 
FACT  =  l.'j 

IF  <  K  .LE.  1  )  PLTURf> 

A 1  =  1.0 
DO  1  I  -  .?*  N 
A 1  -  I 

A  1  =  A 1  *  AT 

CONT  INilf 
FACT  =  A 1 
RETURN 
END 

FUNCTION  AS  INC  *  ) 

ASTN  =  MAM  X  )  S  iVTC  1.0  -  X*X  )  ) 

RETURN 
END 

SUPRUUTINE  CHA<M*(K»  CH«  P  ) 

THIS  S  U  r-  R  0  li  T  1 M  C  COSPUTLS  CHASES’S  A  F  D  t.CUAL  WEIGHTS 

PLK • S  • 


r>  o 


surroutinf  pro^c  k,  w,  p,  ifault  > 


c 


REAL  C‘  <20.20>.P<?0,20> 


l\'Tl-’Ll?r  matrices  t r>  te^o. 


DO  ?  t  ,i  =  1,21 
00  1C  l  =  1  '  S’1 
C  H  <  I  »  J  )  =  0.0 
PCI. J)  -  0. c 
If  CONTINUf 
?  .  )  C  0  N  T  I  l\i  ‘ r 

CH  ( 1  ♦  1 )  =  1.0 
CH<1,2>  =  0.r- 
CMC?. 2)  =  C.  mb 
PClfl)  -  1.0 
P  C 1  *  2  )  =  3.0 
PC2.2)  -  0.r 
DO  40  ^  =2.  K 
J1  -  J  *  1 
v  =  2  *  J  -  1 

Y  =  "  •  J 
U  =  J 

V  =  J  ♦  1 

CHCl ,J1  >  =  V  /  Y  «  CH C 1 »J> 

CHCJi.^u  =  l.o  /  y  *  rh(j.j) 

PC  l*ill  >  =  U  /  V  ♦  PC1,.J> 

PCJl  .J1  )  =  1.0  /  V  *  PC  J,J) 

no  sc  r  =  ?.  j 

IA  =  I  -  1 

CH<I,J1>  =  Cl. 3  /  Y)  *  CHCT  A  *  J  >  ♦  <*/¥>*  CHCI.J) 

P  ( I  .  J 1  )  =  <1.0  /  V)  *  PCIA»J>  ■*  (U  /  V)  *  »CI,J> 

30  CONTINl 
4  0  COM  I  Ml  ‘ 

RETURN 

FNO 

SUBROUTINE  PA  PR  Y  (  K ,  '»'*  PA,  CH,  P  ) 

r* 

C  THIS  S'.'t  h  Cl.  TIME  COMPUTES  THf  APPROXIMATE  PLK*S 

C 

REAL  U  (  2  0  W  F  A  ( 2  0  ,  2  0  )  ,  P  *  (20  *  ?0  >  ,r>PP  ( .>0 , 2  G I  ,  C  IK  20  ,20  )  ,  P<?0 
I  NT f.  SC L  INLPy  (20), A, -3 
ALPHA  -  1.0  /  3.0 

r. 

Z  INITIALIZE  PA,Pb»P«P 


DO  2  0  2  =  l  ,  2  0 
DO  1 3  !  =  1  ,  2  0 
PA(IiJ)  =  3.0 
PRC  I ,J )  -  C  .0 
PRPCT.v!)  =  0.0 

10  CONT INHr 

2  0  CONT  INUr 

C  ntTr-MIAr  HAMMI.IV  AN  1  MINIMUM  Op  -EIGHTS. 

C 

UMAX  =  UC  1  ) 

UMIN  =  -(  1  ) 

00  ?5  I  =  2,  K 

IF  <  W<  I  )  .LT.  WHIN  )  -MIN  =  UC  I  > 

IF  <  'J  (  I  )  .GT.  UMAX  >  UMAX  =  UC  I  ) 

25  CONTINUf 


o  tn 


SUFRCUT  INK  c  R  0 r  F  (  K»  W.  p»  IpAtJLT  ) 


CUT  =  0.6  5  *  W  1 1 N  ♦  0.?f  *  A  X 

r* 

-  DETERM I  'ME  T  t-r  INDICES  OF  THL  WEIGHTS 

r 

DO  U  I  =  1  »  K 

IF  <  W  <  I  >  .LT.  CUT  )  GOTO  5C 
IMnrx<!)  =  i 
GOTO  4’ 

id  INDEX (I >  - 

4?  CONT  l\U  r 

C 

C  COMPUTE  T  HT  NUMBER  OF  LARGE  WT* 

M  r  0 

DO  F  0  I  =  1  .  K 
M  =  m  ♦  INDEX!  T  ) 

G  CONTINUE 

C  IF  ALL  *TS  L  A  R  G '  OR  SMALL  SET  PA-P 
C 

I  F  (<•■*.  r  •  0  )  .AMO  .  ME  .K  )  )  GO  TO  70 
DG  *>f!  L  =  I  *  K 
PA(L,K)  =  F  ( L  t  K  > 

6  f  CONT  I M  D !! 

GO  TO  P70 
C 

C  IF  A  =  0  At  :  P=C  S'.  T  Fr:=PL* 

C 

7T  N  =  IMriFX(l  )  ♦  TNDL’Xm 
IF  <t.  ..ME.  D  GO  TO  90 
DO  8C  L  :  1»  K 
PB<L  »K  )  =  P  (It1'  > 
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